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^ . Abstract 

•/^ ' In this paper we discuss constraints on two-dimensional quantum-mechanical sys- 

^ ■ tems living in domains with boundaries. The constrains result from the requirement 

of hermicity of corresponding Hamiltonians. We construct new two-dimensional fam- 
ilies of formally exactly solvable systems and applying such constraints show that in 
real the systems are quasi-exactly solvable at best. Nevertheless in the context of 

(^ , pseudo-Hermitian Hamiltonians some of the constructed families are exactly solvable. 
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1 Introduction 

It is well known that exactly solvable systems play very important role in quantum theory. 
Unfortunately number of such systems is quite limited. This considerably narrows their 
applications. Such a situation stimulates interest to quasi-exactly solvable systems jH 13 El El- 
In contrast to exactly solvable models in quasi-exactly solvable systems the spectral problem 
can be solved partially. Nevertheless such systems are very interesting. Besides modeling 
physical systems [5J they can be used as an initial point of the perturbation theory or 
to investigate various nonperturbative effects ^. Furthermore, recently in the series of 
papers [7j (see also Refs. [Hlin]) it was revealed a connection between quasi-exactly solvable 
models and supersymmetric systems with polynomial superalgebras [lOj. Also we can hope 
that progress in understanding quasi-exactly solvable quantum-mechanical systems will allow 
to find out methods of constructing quasi-exactly solvable models in quantum field theory. 
The paper has the following structure. In section El a brief introduction into the Lie- 
algebraic approach to constructing quasi-exactly solvable systems is given. In section |3] we 
discuss constraints on wave functions of systems living in domains with boundaries. It is 
shown that for corresponding Hamiltonians to be Hermitian the wave functions must have 
a specific behaviour at boundaries. The role of the constraints is illustrated by an example 
of a known two-dimensional quasi-exactly solvable system. In section |3] we construct new 
families of formally exactly solvable systems. Application of the restriction on the behaviour 
of wave functions at boundaries leads to the conclusion that the systems are quasi-exactly 
solvable even if the wave functions are normalizable. Discussion of results is presented in 
section El 

2 General aspects of quasi-exactly solvable systems 

From the general viewpoint an operator is quasi-exactly solvable if it has a finite-dimensional 
invariant subspace. Bearing in mind applications to the quantum mechanics we are interested 
in second order differential operators, Hamiltonians. The most famous method of construct- 
ing quasi-exactly solvable differential operators is the Lie-algebraic approach P1I21 El El • In 
this section we briefly discuss of main ideas of the approach. 

The one-dimensional and multidimensional cases have to be treated separately. The one- 
dimensional case is the most elaborated one. The simplest finite-dimensional subspaces are 
spaces which admit the following monomial basis: 

Tn = span {1,2;, 2;^...,2;"} . 

For such spaces it is not difficult to find invariant (quasi-exactly solvable) differential oper- 
ators of the first order: 

J_ = — , Jo = AT^ _ !iZLl, J+ = z{N,~n), (1) 

dz 2 

where A^^ = z-j^. These operators form finite-dimensional representations of the algebra 
s[(2,M) for n G N. Any second order operator taken as quadratic combination of the first 
order operators, 

H' = C'jaJt + C'Ja = - P^i^)^ + • • • , (2) 



is automatically quasi-exactly solvable^. Here P^iz) is a polynomial of order 4. In the 
one-dimensional case any second order operator can be reduced to the Schrodinger form: 

/dz (P 

, ^ H= -—^ + V(x), 

where \l/(x) is a wave function of the Hamiltonian in the canonical form. This is essential 
difference of the one-dimensional case from multidimensional systems. 

The complete classification of linear differential operators with invariant subspaces of 
monomials is given in Ref. ||12J (more recent discussions see in Ref. tl3|). 

In multidimensional case the situation is much more complicated [2]. Nevertheless the 
idea is the same. We construct finite-dimensional representation of some Lie algebra in terms 
of first order differential operators [3 E] ? then one can look for Hamiltonians considering 
quadratic combinations of the first order operators: 

H' = C'^'JaJb + CMa = - g^" {V, - A,) (V. - A,) + V, 

with 

V = g^^'A.A, - A^;^, 

where V^ and " ; //" stand for the covariant derivative corresponding to the metric g^^ and 
Such a Hamiltonian can be reduced to the Schrodinger form only if the self-consistency 
conditions are satisfied: 

d^A, - d,A^ = 0. (3) 

From another point of view these equations are the necessary condition for the Hamiltonian 
to be Hermitian. The general solution of the constraints is unknown! Therefore in the 
multidimensional case the general classification of quasi-exactly solvable linear differential 
operators is unknown even for invariant subspaces of monomials. However as it was pointed 
out in Ref. ^ the problem with resolving the constraints can be avoided if one passes 
from D- to {D + l)-dimensional system. Details of the procedure can be found in Ref. P]). 
Nevertheless here we do not adopt this standpoint. 

3 Boundary conditions 

When constructing quasi-exactly solvable operators in the Lie-algebraic approach very often 
resulting systems live in domains with boundaries.^ If so one has to take into account 
behaviour of wave functions at the boundaries in addition to their normalizability. 
Generally the system must evolve in a domain with a positively defined metrics 

5={x'^|(7n>0Udet||^^,|| > O} . (4) 

Here we investigate situation when components of the inverse metric g^^ are polynomial. 
Therefore it is more convenient to define the boundaries of the domain (j3)) as roots of the 
inverse metric determinant: 

dS d {x^'\ detll^'^'^ll = O} . 

^In principle one has to proof that the finite-dimensional subspace corresponds to wave functions with 
finite norms. In the case of the algebra s[(2,M) this question was completely investigated in Ref. |ll| . 
^Here we imply the multidimensional case. 



To provide hermicity of the corresponding Hamiltonian the following conditions have to be 
imphed: 

^g^'ifd^ij = (5) 

dS 

for any f{x) and ifj{x) from the domain of the Hamiltonian. 

In the two-dimensional case the boundary can be locally given by the equation x = ^{y). 
Then from the conditions (0) one can infer that the wave functions from the domain of the 
Hamiltonian have the following behaviour at the boundary: 

'^ ^ g~^ \x - ^iy)r with tt > -. (6) 

The normalizability of such a wave function leads to the inequality a > — | . So the hermicity 
implies the more strict inequality but often it is not taken into account. For example, consider 
the following Hamiltonian [2..: 

-H' = x{l + x)dl + yil + y)dl-2xydl 

+ {l + x){l- ex) d^ + {l + y) (1 - cy) dy + 2c {jx + jy) . (7) 

It can be represented in terms of generators of the algebra 5[(2,M) © s[(2,M). In this case 
the inverse metric is 

u auu ^ fx{l+x) -Xy 

V ~^y 2/(1 + 2/). 

Its determinant has the form 

detllsf^^ll = xy{l + x + y) . 

Therefore we can consider the system in the domain S* = | (x, ?y) I x > U y > O} . In principle 
there are four domains with positively defined metric but the consequence is the same for 
all of them. The wave functions corresponding to the quasi-exactly solvable sector have the 
structure 

where Pol{x, y) is a polynomial in x and y. One can see that such functions are normalizable 
in the domain S for c < 0. However functions of such a form do not belong to the domain of 
the Hamiltonian since they have improper behaviour at the boundaries. Therefore, actually 
the Hamiltonian (|7j) is not quasi-exactly solvable. Nevertheless it can be of interest in the 
context of pseudo-Hermitian Hamiltonians ^H] (also see the discussion below). 

4 New families of quasi-exactly solvable systems 

Formally exactly- solvable systems 

Exactly solvable systems form a subset of quasi-exactly solvable systems, because evidently 
they have (an infinite flag of) finite-dimensional invariant subspaces. 

Let us start discuss of this point in the context of the Lie-algebraic approach from the 
one-dimensional case (Q), 0. One can select two exactly solvable operators of the first order. 



^ and z-^. A Hamiltonian constructed in terms of these operators has the following general 
form: 

The resulting Hamiltonian is exactly-solvable, but only formally. To proof the exact solv- 
ability one has to check normalizability of corresponding wave functions. Not all systems 
given by the Hamiltonian are exactly solvable in this sense. 

Now we pass to the two-dimensional case. The following first order operators 

Ny = ydy, N^ = xd^, Lq =dy, Lp = y^d^, (8) 

where p = 0, 1, . . . , m, are exactly solvable because they are invariant on the spaces 

^„,„= span {xV} (9) 

ma + h^mn 

for any n G N. For any fixed n it is possible introduce yet another operator 

Lni,n = y {mN^ + Ny- mn) , 

which is quasi-exactly solvable. All these operators form a series of representations of the 
algebra R™"^^ ^ g[(2,M) J^. It is worth noting that the case m = 1 has to be considered 
separately since in this case there exist additional quasi-exactly solvable operators of the 
first order, e.g. see Ref il6|. 

The general form of the Hamiltonian constructed in terms of the first order operators (jHI) 
is 

-H= {pox' + Pm{y)x + P2m{y)) dl + [PM^ + P„.+l (?/)) 9^^ + P2(?/)aJ 

+ [qox + Qm{y))dx + Qi{y)dy, 

where Pq, Qq G M, Pq{y) and Qp{y) are polynomials of degree q and p with real coefficients. 
The Hamiltonian ()10|) does not depend on the integer parameter n and, as a consequence, has 
an infinite fiag of finite-dimensional invariant subspaces (jH)) parametrized by n. Therefore, 
formally it is an exactly solvable operator J17J. 

From the expression (|l(Jj) we infer the form of the inverse metric: 

pun _ t Pox^ + Pm{,y)x + P2m{,y) Pi{,y)x + Pm+iiy) 



'^' " ' PM^ + Pm+M P2{y) 

In general the system given by the Hamiltonian (|iup is too complicated to resolve the 
constraints (jH)). For simplicity we investigate the cases of factorizable and unfactorizable 
metric determinants: 

det Wg'-'^W ~ (x - ^y)) {x - ^y)) and det H^^^l ~ (^ - 6(y))' + 6(l/)'. 



To this end we have to fix the coefficient functions, e.g.: 

P2{y) = P2Pi{y)\ 
Pm+i{y) = ^Piiv) {P2PUy) + iPoP2 - 1) iUy) + Uy)) 

P2m{y) = Y (Pm{y) + ^^^^]r^Hy) + ^2(2/))) + ^^^fr^UyMy) 

for the fist case and 

P2iy) = P2Pi{y)\ 
P„+i(y) = Pi{y) {j-Pm{y) + boP2 - 1)6(1/) 

P.n.{y) = P2 (^F.(y) + ^-^^^m) ' + ^^^ Hyf + Uyf) 

for the second case. The functions ^i{y) are polynomials of degrees no more then m. 

By shifting and rescaling the variable y we can fix the form of the linear function Pi{y)- 
There are two different cases. Let us discuss the first one: 

Pi{y) = y- 

In this case there are several solutions of equations (jH)) with factorizable and unfactorizable 
metric determinants. For all of the solutions the coefficient functions Qi{y), Qm{y) and 
Pm{y) are given by 

Qi{y) = Qiy, 

^ , , Pb^CJy) + QiP2yqL{y) - ((2po - go - 2) P2 + 2Qi)g^(y) 

^"^^^^" l-(2po-go-i)P2-gi 

+ (po-go + ^-i)(6(2y) + 6(2/)), 
"^-^'^ = \-{2n-q.-i)P2-Q. + — f^(^^^^) + ^'^'^^- 

According to the first solution the polynomials ^i{y) are arbitrary, 



(11) 



1-g 

and the wave functions have the form 



qm{y) = ( p - 2po + qo + 1] ^i{y) 



^ = g-^\ynx-Uy)fPol{x,y), (12) 

where Pol{x, y) G J-'m,n i = 1 or 2 and 

^ QiVq + Po - go - 1 ^ goA + P2 - Qi - 1 

"" 2(poF2-l) ' ^" 2(poP2-l) 




Uy) 



^liy) 




6/" + fi(y) 



fi(>') 




Figure 1: Schematic representations of boundaries. 



The determinant of the inverse metric has the factorizable form 



det 



1^' 



l^u I 



y'{^-Uy)){x-Uy))- 



Possible configuration of boundaries in this system is schematically represented on the left 
plot of figure HJ From the structure of the wave functions (fT^ one can conclude that for the 
Hamiltonian to be Hermitian the system should live in the domain II (//?) with x ^ ^2(2/) 
and y ^ {y ^ 0) or IIl (IIr) with x ^ Ciiv) and y ^ {y ^ 0). But it is easy to check 
that for any choice of the parameters a > | and (3 > \ the wave functions (fT^ are not 
normalizable. Therefore the system is actually not quasi-exactly solvable. 

In the regions IIIl and IIIr the wave functions (fT^ do not belong to the domain of 
the Hamiltonian. To consider these domains it is necessary to discard the hermicity of 
the Hamiltonian. Nevertheless in the quasi-exactly solvable sector matrix elements are real 
hence corresponding eigenvalues are real or form complex conjugated pairs. But this pattern 
is specific to pseudo- Hermitian Hamiltonians jT^I. Therefore the system in the domains IIIl 
and IIIr can be interesting from this point of view. Moreover some choice of the parameter a 
can provide power-like falloff and normalizability of the wave functions ()12|) (at least part of 
them). In this sense the system is quasi-exactly solvable. Besides if we discard the hermicity 
then it is possible to consider another configurations of boundaries when the corresponding 
system is quasi-exactly solvable or even exactly solvable. Such a configuration is represented 
on the right plot of the figure (Q). Indeed in this case the system living the finite region / 
is exactly solvable because wave functions are normalizable for any values of the parameters 
a > — I and P > — |. 

According to the second solution 



6(?/) = ei(?/)+6/ 



Qmiv) 



KlV" ~ ^^° + "^^ + V ^'^^^ ^ "^"^^ 



and the wave functions have the following structure 

'0 = 9~^ \yr\x-Ciiy)f \x-^iiy) -i2y^\' Pol{x,y), 

where Pol{x,y) G J-'m^n, k = Pj^^ ^ N and 



(13) 



QiPo + Po - ?o - 1 , a {Q0P2 + P2-Q1- 1)6 - P2qv 

Ot = — ^^ 7^ 1, p = 



2{poP2 - 1) 



2{poP2 - 1)6 



1, 7 



P2qn 



2(^0^2-1)6' 



The determinant of the inverse metric is factorizable 

det Wg^^W ^y^{x- ^y)) (^ " Uv) - 62/') • 

Possible configuration of boundaries in this system is schematically represented on the central 
plot of the figure ^ From the structure of the wave functions (fT!^ one can conclude that 
for some choice of the parameters the wave functions are normalizable with power-like falloff 
and the system is quasi-exactly solvable in any of the domains. 

According to the third solution the polynomials iiiy) and Qmiy) are not fixed, Qi = 
1 — P2 (2po — Qo — ^) and the wave functions have the following form 

X go — 1 

V' = 9^"^ \y\~ ^° Pol{x, y), 

where Pol{x,y) G J^m,n- Such functions do not belong to the domain of the Hamiltonian or 
are not normalizable for any choice of the parameters. Therefore the corresponding system 
is not quasi-exactly solvable. Nevertheless if one discards the hermicity of the Hamiltonian 
then for the situation represented on the right plot of figure ^ the system can be quasi- 
exactly solvable (in the regions IIl and IIpi) or even exactly solvable (in the region J) for 
some choice of the parameters go and Pq. 

Now we pass to the metric with the unfactorizable determinant, i.e. it has no roots in 
the variable x. In this case there is the following nontrivial solution: 

6(2/) = ^i{y) + 6?/^ qm{y) = qmi2y^ - {Qik -k + 2po-qo-i) 6(2/) 

and the wave functions have the following structure 

^ = r ^ \yr ((^ - Uy)? + e^y^'Y e^^-*^^^Fo/(x, y), (14) 

where Pol{x, y) E J-'m,n, k = P2^ E N and 

_ k{Qipo+po -go - 1) _ Qik + k - gp - I 1 _ Qm 

"" 2ipo-k) ' ^" 4(/c-j9o) 2' ^~2{k-poy 

The determinant of inverse metric has the form 

det 11^^11- 1/2 ((x- 6(2/))' + e2V'). 

The determinant vanishes only at the line y = 0, therefore the system can live on the left 
or right semiplane of Oyx. One can see that negative values of the parameter (] with large 
enough modulus can provide power-like falloff and normalizability of a finite number of the 
wave functions (|T4J1 . Therefore the system is quasi-exactly solvable. 

Now let us consider the second possibility of fixing the linear function Pi{y): 

Pi{y) = Pi- 



In this case for all of the solutions the coefficient functions Qi{y), Qmijj) and Pmijj) are given 
by 

^ I , _ P^Ph'LJy) + PiPiQiq'mJy) - ^2 (2po - go + 2Q1) g^y) 
^"^^^^" l-(2po-go + gi)P2 

+ (Po-go + Qi)(ei(2/) + 6(?/)), 

1 - (2po - go + <5i) P2 P2 

According to the first solution the polynomials ii{y) are arbitrary, 

(lm{y) = f -^ - 2po + go - Qi j 6(1/) 

and the wave functions have the form 

i^ = g-'*\x~^,{y)re-('yPol{x,y), (16) 

where Pol{x,y) G J-'m,n, i = 1 or 2 and 

P2 (go - Qi) - 1 . . go-Po(i'2Qi + i) 

a = — — — — 1, p — 



2ipoP2-l) " 2P,{poP,-l) 

The determinant of the inverse metric has the form 

det\\g''^\\r^{x-Uy)){x-^2{y)). (17) 

For the Hamiltonian to be Hermitian the functions ^i{y) have to obey the inequality C,i{y) ^ 
^2(?/) or Ci{y) ^ C2{y)- Such a configuration is represented on the left side of figure [T] One 
can check that for any choice of the parameters a and f3 the wave functions (fT^ are not 
normalizable. Therefore the system is actually not quasi-exactly solvable. Nevertheless if 
non-Hermitian Hamiltonians are allowed then for configurations of boundaries represented 
on the right plot of figure^ the systems are exactly solvable in the regions /, IIl and IIr. 
Also there is a solution with the determinant of the form ()17|) with arbitrary polynomials 
6(?/); Qmiy) and wave functions with the structure ip ~ g~^/'^e~^'yPol{x,y). In this case 
hermicity of the Hamiltonian or normalizability cannot be provided. Therefore the system 
is not quasi-exactly solvable. But if we discard the hermicity of the Hamiltonian then for 
the configuration represented on the right plot of figure [T] the system is exactly solvable if it 
lives in the domains /, IIl and IIr. 

5 Conclusion 

In this paper we have considered quantum-mechanical systems in domains with boundaries. 
Such a situation is usual for the Lie-algebraic approach fll El to construction of multidimen- 
sional quasi-exactly solvable Hamiltonians. The hermicity of the operators prescribes the 

9 



specific behaviour of wave functions which belong to domains of Hamiltonians. By an exam- 
ple it was shown that not all of known quasi-exactly solvable systems are in fact quasi-exactly 
solvable even if corresponding wave functions are normalizable. 

Besides we have constructed new two-dimensional families of formally exactly solvable 
systems. Application of the restrictions the behaviour on wave functions in a domain with 
boundaries leads to the conclusion that the constructed systems are quasi-exactly solvable 
at best. 

If one discards the hermicity then the behaviour of the wave functions is governed only by 
the normalizability. It is worth noting that the corresponding matrix elements are real. In the 
quasi-exactly solvable sector the corresponding matrix is finite dimensional therefore quasi- 
exactly solvable eigenvalues are defined by an algebraic equation with real coefficients. This 
means that the eigenvalues are represented by real values and complex conjugated pairs of 
numbers but this pattern corresponds exactly to the spectrum of a pseudo-Hermitian Hamil- 
tonian ^3]. Therefore the formally quasi-exactly solvable and exactly solvable systems can 
be interesting in this context. Of course, the exact relationship between such quasi-exactly 
solvable operators and pseudo-Hermitian Hamiltonians requires more detailed investigation. 
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